Disordered XY models and Coulomb gases: renormalization via traveling waves 
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We present a novel RG approach to 2D random XY models using direct and replicated Coulomb 
gas methods. By including fusion of environments (charge fusion in the replicated CG) it follows the 
distribution of local disorder, found to obey a Kolmogorov non linear equation (KPP) with traveling 
wave solutions. At low T and weak disorder it yields a glassy XY phase with broad distributions 
and precise connections to Derrida's GREM. Finding marginal operators at the disorder-induced 
transition is related to the front velocity selection problem in KPP equations yielding new critical 
behaviour. The method is applied to critical random Dirac problems. 



Two dimensional random systems have attracted con- 
siderable recent interest in domains ranging from local- 
ization in quantum Hall systems to vortices in supercon- 
ductors. In the context of localization, progress was made 
to characterize the multifractal statistics of 2D wavefunc- 
tions using random Dirac models ||l|,||, extending previ- 
ous studies in ID On the other hand the glassy prop- 
erties of vortex phases with disorder was investigated us- 
ing random XY models. While these lines of studies have 
developed in an apparently disconnected way, they led to 
similar proposals that remarkable connections ex- 
ist between the large fluctuation properties of these sys- 
tems and Derrida's random energy models GREM 
To study these connections further, consistent RG tech- 
niques are needed. Our aim is to develop such an ap- 
proach, which, as in other glassy systems, e.g. in ID 
requires a proper treatment of broad distributions. 

Here we focus on random gauge XY models and discuss 
at the end related random Dirac problems. Recently the 
phase diagram predicted long ago in |^ was reexamined 
using energy arguments. It was proposed that de- 
fects (vortices) are not thermalized at low T, leading to 
a non-reentrant transition line between the XY and the 
defective phases, and to a failure of the conventional 
Coulomb Gas perturbative expansion of This mo- 
tivated several interesting proposals for RG procedures 
P, p]JT^ . However, these approaches, while giving the 
correct topology of the phase diagram, are not fully con- 
sistent, since they do not take into account the renor- 
malization of local disorder. As shown here, this changes 
quantitatively the results of |^,|ri|,|l2| for the XY phase 
and becomes crucial at the disorder driven transition. 

In this Letter we reexamine the RG procedure for the 
disordered Coulomb gas (CG) and for the random gauge 
XY model and propose an approach which allows to treat 
the probability distribution of the local disorder. Techni- 
cally this amounts to introduce composite charge fugac- 
ities and to study their fusion (in the CG sense) upon 
coarse graining, while only dipole fugacities were consid- 
ered previously P, |rT| , p^ . A precise connection between 
the fugacity distribution and the GREM free energy dis- 
tribution is found via the Kolmogorov (KPP) equation 
p3[ which arises as its RG flow equation. Universality 
in the corresponding non linear front selection problem 



unexpectedly translates into the RG universality around 
the disorder driven transition. This allows to describe 
the phase dominated - and the transition driven - by 
rare configurations of frozen defects, where correlations 
are broadly distributed. Restriction to the single charge 
sector yields a RG derivation of the multifractal proper- 
ties of the critical Dirac wavefunction. 

The 2D square lattice XY model with random 
phases M is defined by its partition sum Z[A] ~ 
U^ d0, e-^«[»'^l with: 
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13 J, P = l/T. The 



Aij are independent gaussian random gauge fields, with 



Af- — ncr. This model can be transformed exactly [ni 15 



into a CG with integer charges defined on the sites r of 



the dual lattice with Z[V] = i ^-^^ and 
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-K ^ rirGr 
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where Gj^^ = i(2 — cos(fca;a) — cos{kya)) is the lattice 
Laplacian. The bare disorder potential, 14 = -^Grr'(V x 
A)ri is gaussian with logarithmic long range correlations 
Vkl^-k — 2aK^Gi^. The usual continuum approxima- 
tion with (integer) charges of hard core a and fugacities 
y = e'^''^ of this lattice model is obtained by using the 

asymptotic form Grr' ~ 0-'n\r — r'\/a + j){l — 5rr')- Here, 
the perturbative expansion of Z[V] in y, valid in the di- 
lute limit uniformly over the system, fails. 

Before turning to the more systematic replica ap- 
proach, let us first sketch the direct RG method suited 
to the present case where disorder favors some regions, 
resulting in a site dependent local fugacity y^. Our ex- 
pansion captures the limit where the fugacity is negli- 
gible almost everywhere except in a few rare favorable 
regions. This is achieved by following the local disor- 
der distribution which is not gaussian, a novel feature 
from all previous approaches. We find that the disor- 



der Vr 



naturally splits into two parts, a 



long range correlated gaussian part V/' with logarithmic 
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correlator {V> — V^)^ — AaK^ln{\r — r'\/a) and a lo- one gets 



cal non gaussian part Vr which defines the local fugacity 
variables Zj. — yrexp(±Vr) for ±1 charges which 
have only short range correlations. The RG equation 
for the distribution P(z+,z-) of local environments is 
obtained from two contributions (i) "rescaling": upon 

produces a gaus- 



di 



coarse graining a — > a = ae' 
sian additive contribution to v 



from (V> - V>) 



>\2 - 



) + dl] = {V> -V>f +idv, 



dV-r 



,v Kdl±dv, with dVrdVr' = 2aK^dl5r 



(ii) "fusion of charges" (fusion of environments) upon the 
change of cutoff, two regions with fugacities Zj- , z^ are 
replaced by a single region at f = ^(r' + r") of effective 
fugacities z± — [zl^ + +zLz!j_ ) obtained 

from the relative weight W+/Wo of a charge 1 configu- 
ration (either in r' or r") versus a neutral one (either no 
charge or a dipole). (i) and (ii) yield: 



diP{z+, z_) =OP- 2P(z+, Z-) + 2 ( 6{z+ 



z'+ + z'_l 



1 + zLz'i + . 



z_ + z': 



1 + z'_z'_ 



+ z'^z'L 



P'P" 



(3) 



where {A)p,pn denotes A P{z'_^, z'_)P{zl, z'_L) 

and O = K{2 + z+d,^ + z^d;,_) + aK^{z+d,^ ~ z^d.y 
is the diffusion operator p^ . 

To put this derivation on a firmer footing and to cap- 
ture broad distributions of local fugacities we introduce 
i pTf an expansion of physical quantities in the number of 
points ||l8|, which for the free energy F[V] ~ —TlnZ[V] 
as a functional of the disorder reads: 



\V] 



(4) 



where by definition /rV, .r^ depends only ||T^ on V{r) at 
points Ti, i = l,..fc. The continuum limit of (^) con- 
sists in replacing the sums by integrals with hard core 
constraints around each r,-. The first two terms read: 



-l3f^3=Hl + Wr.r'), 



Bf^^^ =ln( — 



Wr.r' — Wr,r' 
1 + Wr'^r" + VKr,r' + 



Wr 



(1 + Wr,r'){l + Wr',r")(l + W^r.r") ' 



(5) 
(6) 
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z,z_ e 



where lUr.r 

the Boltzmann weight of a dipole. The first term (||) cor- 
responds to the independent dipole approximation | [T^ ] 
while (^ takes into account contributions from triplets 
of sites. Though there are no actual configuration with 
three charges in a given environment (neutrality), this 
term (and higher orders), overlooked in previous ap- 
proaches, is crucial for the renormalization as it leads 
to fusion of environments when coarse graining. Upon 
increase of the cutoff a — > ae'^' on the continuum version 
of (|^), the k point integral gives a correction to order 

dl to the fc - 1 point integral, e.g. Ja<\r'-r"\<ae-" /r,r',r" 

corrects fl^ (with f = i(r' + r")) as 5ln{l + VFr,f) = 
d/[ln(l + iyr,f)-ln(l + T4^;_f)-ln(l + T4^;;f)]. First VFr-r" = 
z'+z'f_ + z'^z'f_ since V^^ = V^,+0{dl) (we denote z^ = z± 
etc.), thus the combination {wry + Wry)/{1 + Wr'r") 
in (H) can be rewritten as lir,} = z+z-c^^^ -Gr,? 
terms of the new fugacity z± defined above. Second, 
f = z+z'_e~^'- """^f -Gr,f (similarly ^ using z^). 



Averaging this correction over disorder yields a result 
corresponding to the RG equation for P (||). A com- 
plete derivation of (H) involves a similar procedure for all 
moments of F[V] |0, provided one adds to (||) the free 
energy sum of all degrees of freedom eliminated up to 
scale I with diFo = -T(ln(l + z'_^z'f_ + z'_zl))) p>p,, . 

Finally, the renormalization of K and a is 
obtained from screening [^|JTl|,[l^: diK~^ — 
-2TT'^{nonp=a)c,dia = -2T:^{no){np^a)- Expanding in 
the number of sites using (i), we obtain [0 the follow- 
ing, which together with (5) forms our complete set of 
RG equations [HPI 



dK- 



dl 



Ad' 



z'.z'L + z'_z'l 



- Az'^z"_z'_z'l_ 



[l + z'+z'L + z'_zlf 



P'P" 



da 



{z'+z'f - z'_z\f 

{I + z'+z'f + z'_z'iy 
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The combinatorics necessary to this method is much 
easier performed using replicas. We start again from (|^) 
and represent as the partition sum of a CG with m- 
vector charges living on the dual lattice sites. Averag- 
ing over disorder, and taking the continuum approxima- 
tion we obtain the m-vector (hard core) CG of partition 
sum expanded in power of the vector fugacity Yn'. 



Z'" = 1 



E 



E 



• -^n 



p>2 Hi. .Up ^1 



n 



with Kbc = KStc — a-K^, all integrals being restricted to 
|ri — > a, and the sum is over all distinct neutral con- 
figurations J2t = 0- is a function of the replicated 
charge n — {n^ , ..n™) with bare value y„ ^--tnbK'"'na ^ 
Since Kb^c 7^ 0, one cannot restrict to single non zero 
component charges as it leads to the erroneous re- 
sults of at low temperature. However, we stress that 
this quadratic form for Y^, which results from the Gaus- 
sian nature of the bare disorder, is not preserved by the 
RG as shown below. We now perform the RG analysis of 
the m-vector CG, extending the scalar case leaving 
the above form unchanged with p3| : 
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(8) 



Rescaling and annihilation of opposite replica charges 



separated by a < jr^ — rj| < 



gives the first term 



of (g) and (g). The second term of @ which comes from 
fusion of two replica charges as usual in vector CG, was 
absent in but is necessary for consistency of RG 

to order Y^. 

Why should one consider the expansion in Yn ? Tech- 
nically, it is valid, together with (pp|), in the limit of 
a small density of vector charges pOf, which here corre- 
sponds to a small density of favorable local regions. In- 
deed we checked that this expansion is identical term by 
term, for m ^ 0, to the expansion in number of points 
of the free energy (||). Thus the set of Yn should en- 
code the full scale dependent distribution P{z±) of local 
disorder, the perturbative parameter being P{z^ ^ 1). 
Remarkably, the correspondence between z^) and 

Yn emerges when performing the analytical continuation 
m ^ of (HH) which we now present. To capture the 
most relevant operators it is sufficient to consider Y^ with 
= 0, ±1 in each replica which, using replica per- 
mutation symmetry, depends only on the numbers n± of 
±1 components of n. This leads to the general parametri- 
sation in term of a function z-): 



where (..)$ = /^^ After some combinatorics 

pl| the limit to — > of can be rewritten equivalently 
as an equation for <I>, detailed in The first term 

in gives a diffusion contribution (2 + O)'^ and the 
second term in yields a term of fusion of environ- 
ments, analogous to the one in (^. From this equation 
TV = ^^^p $(z±) is found to satisfy diN = 27V-c^A/'^ 
and thus converges quickly towards M* = 2/d. We 
thus define the normalized P = which satisfies (^) 

and is naturally interpreted as a probability distribution 
p^ . Finally, with the same definitions, (^) yields (0). 



^(z) -X,=E,-m, 




Front region- {T 

FIG. 1. Scale dependent distribution Pi{z) and its two tail 
regions for T < Tg. Inset: phase diagram 



We first study numerically the RG equations (0J|) and 
find at low T, a < ac an XY phase as in Fig.l {K, a 
converge to K^, an). The typical z goes to zero but P 
develops a broad tail up to ~ 0(1). While in this phase 
and at criticality the concentration of rare favorable re- 
gions decreases, it eventually increases at large I in 
the disordered phase Uc ^ tr. 

To go beyond numerics, we argue |25[| that it is con- 
sistent to discard the z'_^z'!_ terms in the denominators 
in (H). This leads to a closed equation for a single fu- 
gacity distribution P{z) — P{z, z-) which, using the 

parametrization Gi{x) = 1 — (exp(— ze"'''-^"^'-'))^ 
where Ei ^ £Jidl, (3^ 1/T, can be rewritten as: 



-diG ^ DidlG + {I - G)G 



(10) 



The diffusion coefficient is Di = \uiJf and by construc- 
tion G/(— oo) — 1 and G'/(-|-oo) — 0. Remarkably, for con- 
stant D this is the much studied KPP equation, which de- 
scribes diffusive invasion of an unstable state (G — 0) by 
a stable one (G = 1), also related to branching diffusions 
and glassy REM- like models [|6). It is known ||l3|,|^ that 
Gi{x) converges at large I towards traveling waves solu- 
tions h(x — mi) selected by the behaviour at infinity of 



G/=o(^) 



This implies Pi{z) — > z p(lnz — f3Xi) 



with Xi— mi — El {Xi < in the XY phase, see Fig.l). 

At low T in the XY phase, Pi{z) becomes very broad 
and one must distinguish two different tails. As shown 
in Fig.l, the bulk of the distribution (typical values) 



is located around ztyp 



It corresponds to the 



y ' 

SDl ■ 



i+A") within the tail of 
Thus for T < Tg the 



front region which has a tail of size ^/l ahead of the 
front. There from the velocity selection studies ^3 27 
for T > Tg = Jy^ a/2 we find the front position mi ~ 
2{(3^^ +D(3)l. For T < Tg the velocity freezes with to/ = 
^/D{Al-^\Ill + 0{^)) and hiiy) - ^ U. 

This corresponds to Pi{z) ^ z^' 
the front, with fj. = T/Tg < 1. 
distribution function of In z travels at the relative ve- 
locity diXi — J{\/8a — 1), which determines the phase 
diagram: it is negative (decrease of P(l)) in the low T 
XY regime, positive for a > cFc- Furthermore, at low T, 
there is also a far tail ahead ~ / of the front which cor- 
responds to rare events z ~ 1, of small probability -Pz(l), 
but which dominate average correlations (and thus diK 
and dia). The linearized KPP equation, valid in this re- 
gion, leads to Pi{z) - P/(l)z-^+'^ with Ji = T/T* < 1 
and to Pi(l) - e(2-i)', for T < T* = 2aJ. A more 
detailed study of the XY phase is given in |l^ . 

In the high T regime of the XY phase, Pi{z) is not 
so broad, and one recovers from ( p^ ) the usual RG 
result diy — (2 — K + (jK'^)y for the average fu- 
gacity yi = {z)i < +oo (~ ztyp for T > Tg), using 



Giix) 



-P(x-Ei). 



z) p, at large x. 
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The critical behaviour at the transition from the XY 
to the disordered phases is determined by the front re- 
gion, since the velocity diXi vanishes. Here we sketch 
the analysis at T = 0: defining z — e^^ and using 
lim^^oo fipi{l3{u - Xi)) = -h\{u - Xi), (0) yields: 



diJ- 



8d' 
IF 



h[{^2Xi) ■ dia 



IF 



H-2Xi 



(11) 



(which is of order -P(l)^ perturbative for Xi ^ 1). 
From the universal corrections [^3[^ to the velocity: 
diXi = J{V8^ - 1) - 3VD/21 + 0(ri) we obtain a 
projection of the RG flow on the plane tr ~ Uc = | and 
g ^ ^ h{-Xi) which reads p|: 



dig = (16(cr- (Tc) - 



di{a - (Jc) = g^ 



(12) 



yielding gi 

1 



e 1° 



I 2 at criticality and a correlation length 
T. This new universality class U% is different 



from KT and from the prediction of jl^,^ • Note that al- 
though most details of Piz±) e.g. its bulk, depend on the 
cutoff procedure (and fusion rule...), here the universality 
appears in a remarkable way. It comes from the indepen- 
dence of the velocity and the front tail (which also deter- 
mine the relevant operators) on the precise form F[G] of 
the non linear term in (|o|) (see Q). 

Finally, our RG also applies to the problem of a single 
charge Z = e~^^ in a random potential with logarith- 
mic correlations, related to diffusion in random media j^] 
or wavefunction of 2D Dirac fermions in a random mag- 
netic field @-||. The same decomposition of disorder, 
and fusion of environments (z = z' + 2", fugacities being 
local partition sums) yields (p^). We recover for Pi{z) the 
mapping to directed polymers (DP) on Cayley trees, con- 
jectured in with the same universal intensive free 
energy and wavefunction multifractal spectrum [p|,pTt. 

To conclude, we developed a RG approach to random 
XY models, disordered CG and random Dirac problems. 
By following the whole fugacity distribution, it appears 
perturbative in the concentration of rare favorable re- 
gions, which corresponds to the vector fugacity in the 
replica method. This expansion is highly non perturba- 
tive in the original fugacity y. A precise connection to 
the free energy distribution of DP on Cayley trees and 
GREM arises from the RG and turns out to be cru- 
cial to describe the disorder driven transition. 

We thank B. Derrida and V. Hakim for useful discus- 
sions about the KPP equation. 
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